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We present details about an experimental method based on 
the Landau Zener model which allows to measure very small 
tunnel splittings A in molecular clusters Fes. The measure- 
ments are performed with an array of micro-SQUIDs. The 
observed oscillations of A as a function of the magnetic field 
applied along the hard anisotropy axis are explained in terms 
of topological quantum interference of two tunnel paths of op- 
posite windings. Transitions between M = — S and (S — n), 
with n even or odd, revealed a parity (symmetry) effect which 
is analogous to the suppression of tunneling predicted for half 
integer spins. This observation is the first direct evidence of 
the topological part of the quantum spin phase (Berry phase) 
in a magnetic system. The influence of intermolecular dipole 
interactions on the measured tunnel splittings A are shown. 



I. INTRODUCTION 

Magnetic molecular clusters are the final point in the 
series of smaller and smaller magnets from bulk matter 
to atoms. Up to now, they have been the most promis- 
ing candidates to observe quantum phenomena since they 
have a well defined structure with well characterized spin 
ground state and magnetic anisotropy. These molecules 
are regularly assembled in large crystals where often all 
molecules have the same orientation. Hence, macroscopic 
measurements can give direct access to single molecule 
properties. The most prominent examples are a dode- 
canuclear mixed-valence manganese-oxo cluster with ac- 
etate ligands, Mni2 acetateEl, and an octanuclear iron(III) 
oxop hydroxo cluster of formula [Fe802(OH)i 2 (tacn) 6 ] 8+ , 
FesB, where tacn is a macrocyclic ligand. Both systems 
have a spin ground state of S = 10, and an Ising-type 
magneto-crystalline anisotropy, which stabilizes the spin 
states with M — ±10 and generates an energy barrier for 
the reversal of the magnetization of about 67 K for Mni2 
acetate and 25 K for Fes- 

Fe§ is particular interesting because its magnetic re- 
laxation time becomes temperature independent below 
0.36 K showing for the first time that a pure tunneling 
mechanism between the only populated M = ±10 states 
is responsible for the relaxation of the magnetizationB. 
Measurements of the tunnel splitting A as a function of 



a field applied in direction of the hard anisotropy axis 
showed oscillations of A, i.e. oscillations of the tunnel 
ratea. In a semi-classical descriptionBQ, these oscilla- 
tions are due to constructive or destructiveAnterference 
of quantum spin phases of two tunnel pathsCHj. Further- 
more, parity effects were observed when comparing the 
transitions between different energy levels of the systemcl 
which are analogous to the parity effeet between systems 
with half integer or integer spinsBLj. An alternative 
explication in terms of intermediate spin was also pre- 
sentedEJ. Hence, molecular chemistry had a large impact 
in the research of quantum tunneling of magnetization at 
molecular scales. 
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FIG. 1. Schematic representation of our magnetometer 
which is an array of micro- SQUIDs. Its high sensitivity al- 
lows us to study single crystals of the order of 10 to 500 /im 
which are placed directly on the array. 

II. MICRO-SQUID TECHNIQUE 

The technique of micro-SQUIDs is very similar to the 
traditional SQUID technique. The main difference is that 
the pick-up coil is replaced by a direct coupling of the 
sample with the SQUID loop (Fig. [lj). When a small sam- 
ple is directly placed on the wire of the SQUID loop, the 
sensitivity of the micro-SQUID technique is tpi orders of 
magnitude better than a traditional SQUIDlI 2 ] reaching 
10~ 17 emu. This sensitivity is smaller when the sample 
is much bigger than the micro-SQUID. 

Our new magnetometer is a chip with an array of 
micro-SQUIDs. The sample is placed on top of the chip 
so that some SQUIDs are directly under the sample, 
some SQUIDs are at the border of the sample and some 
SQUIDs are beside the sample (Fig. |l|). When a SQUID 
is very close to the sample, it is sensing locally the mag- 
netization reversal whereas when the SQUID is far, it is 
integrating over a bigger sample volume. 

The high sensitivity of this magnetometer allows us to 
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study single Fes crystalalil of the order of 10 to 500 /im. 
The magnetometer works in the temperature range be- 
tween 0.035 and 6 K and in fields up to 1.4 T with sweep- 
ing rates as high as 1 T/s, and a field stability better than 
a microtesla. The time resolution is about 1 ms allow- 
ing short-time measurements. The field can be applied 
in any direction of the micro-SQUID plane with a pre- 
cision much better than 0.1° by separately driving three 
orthogonal coilalj. In order to ensure a good thermali- 
sation, the crystal is fixed by using a mixture of araldite 
and silver powder. 
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FIG. 2. Temperature (a) and field sweeping ratepft)) de- 
pendence of hysteresis loops of Fes molecular clusteraU. Res- 
onant tunneling is evidenced by equally separated steps of 
AH Z » 0.22 T which, at T < 360 mK, correspond to tunnel 
transitions from the state M = — 10 to M — 10 — n, with 
n — 0,1, 2.... The resonance widths of about 0.05 T aee.|djie 
to mainly dipolar fields between the molecular clusterstJlij. 

Typical measurements of magnetic hysteresis curves 
for a crystal of molecular Fes clusters are displayed in 
Fig. U and ||[ The field was swept in direction of the 
easy axis of magnetization evidencing about equally sep- 
arated steps at H z w n x 0.22 T (n = 1, 2, 3 ...) which 
are due to a faster relaxation of magnetization at par- 
ticular field values. The step heights (i.e. the relaxation 
rates) change when a constant transverse field is applied. 
It is the purpose of these article to present a detailed 
study of this behavior which is interpreted in terms of 



resonant tunneling between discrete energy levels of the 
spin Hamiltonian S = 10 of Fes- 
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FIG. 3. Hysteresis loops measured in the presence of a con- 
stant transverse field, at 0.04 K. Insets: enlargement around 
the field H — 0. Notice that the sweeping rate is ten times 
slower for the measurements in the insets than that of the 
main figures. 



III. LANDAU ZENER METHOD 

The simplest model describing the spin system, of Fes 
molecular clusters has the following Hamiltoniancl: 

H = -DS 2 Z +E(S 2 X - Si) +gfi B SH (1) 

S x , S y , and S z are the three components of the spin oper- 
ator, D and E are the anisotropy constants, and the last 
term of the Hamiltonian describes the Zeeman energy 
associated with an applied field H. This Hamiltonian 
defines a hard, medium, and easy axes of magnetization 
in x, y and z direction, respectively (Fig. ^). It has an 
energy level spectrum with (25+1) = 21 values which, in 
first approximation, can be labeled by the quantum num- 
bers M = —10, —9, ...10. The energy spectrum, shown 
in Fig. ^, can be obtained by using standard diagonalisa- 
tion techniques of the [21 x 21] matrix describing the spin 
Hamiltonian S = 10. In the low temperature limit (T < 
0.36 K) only the two lowest energy levels with M = ±10 
are occupied. The avoided level crossing around H z = 
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is due to transverse terms containing S x or S y spin 
operators (see inset of Fig. |7|). The spin S is 'in res- 
onance' between two states when the local longitudinal 
field is close to the avoided level crossing (< 10~ 8 T for 
the avoided level crossing around H z = 0). The energy 
gap, the so-called tunnel spitting A, can be tuned by an 
applied field in the rry-plane (Fig. ||) via the S X H X and 
SyHy Zecman terms. It turns out that a field in H x direc- 
tion (hard anisotropy direction) can periodically change 
the tunnel spitting A. In a semi-classical description, 
these oscillations are due to constructive or destructive 
interference of quantum spin phases of two tunnel paths 
(Fig. ^). The period of oscillation is given byO: 



AH = — V / 2E(E + D) (2) 

The most direct way of measuring thc-tunnel splitting A 
is by using the Landau- Zener modelOEa which gives the 
tunneling probability P when sweeping the longitudinal 
field H z at a constant rate over the avoided energy level 
crossing (see inset of Fig. ||): 

77 A 2 , 

M , M / 

PM M> = 1 — e 2ngn B \M-M'\dH/dt (3) 

Here, M and M' are the quantum numbers of the avoided 
level crossing, dH/dt is the constant field sweeping rates, 
g k, 2, [is the Bohr magneton, and h is Planck's con- 
stant. In the following, we drop the index M and M' . 
For very small tunneling probabilities P, we did multiple 
sweeps of the resonance transition. The relaxed magne- 
tization after N sweeps is given by (n = 0): 

M(N) ~ exp [-2PN] = exp [-Ft] (4) 

where N — -jr^jri is the number of sweeps over the 

level crossing, T = 2P\^- = 4r"3W" is tne overa11 
Landau-Zener transition rate, AM is the change of mag- 
netization after one sweep, and A is the amplitude of the 
ac- fieldMl. We have therefore a simple tool to obtain the 
tunnel splitting by measuring P, or M(N) for P << 1. 

In order to apply the Landau-Zener formula (Eq. |^), 
we first saturated the-,sample in a field of H z = -1.4 T, 
yielding M\ n = — M S E3. Then, we swept the applied field 
at a constant rate over one of the resonance transitions 
and measured the fraction of molecules which reversed 
their spin. This procedure yields the tunneling rate P 
and thus the tunnel splitting A (Eq. ||). We first checked 
the predicted Landau-Zener sweeping field dependence 
of the tunneling rate. This can be done by plotting the 
relaxation of magnetization as a function of t = N d ^ dt ■ 
The Landau- Zener model predicts that all measurements 
should fall on one line which was indeed the case for 
sweeping rates between 1 and 0.001 T/s (fig. |). The 
deviations at lower sweeping .rates, are mainly due to 
the 'hole-digging mechanism'Ej which slows down the 
relaxationEil. In the ideal case, we should find an ex- 
ponential curve (Eq. ||). However, this might be only the 



case in the long time regime (see inset of Fig. The 
origin is not clear but it might be related to dipolar inter- 
actions and hyperfine couplings. For comparison, there 
is also a relaxation curve without ac-field at H = which 
shows much slow relaxationta. 
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FIG. 4. Unit sphere showing degenerate minima A and B 
which are joined by two tunnel paths (heavy lines). The hard, 
medium, and easy axes are taken in x, y and z direction, 
respectively. The transverse field H tr ans is applied in the xy 
plane at an azimuth angle ip. At zero applied field, the giant 
spin reversal results from the interference of two quantum 
spin paths of opposite direction in the easy anisotpppy plane 
yz. By using Stokes theorem it has been shownEI that the 
path integrals can be converted in an area integral, giving 
that destructive interference, that is a quench of the tunneling 
rate, occurs whenever the shaded area is kn/S, where k is an 
odd integer. 
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FIG. 5. Zeeman diagram of the 21 levels of the S — 10 
manifold of Fes as a function of the field applied along the easy 
axis (jjj). Form the bottom to the top, the levels are labeled 
with quantum numbers M — ±10, ±9, ...0. The levels cross at 
fields given by fioH n ~ n 0.22T, with n — 1,2, 3.... The inset 
displays the detail at a level crossing where the transverse 
terms (terms containing S x or/and S y spin operators) turn 
the crossing into an avoided crossing. The higher the gap A, 
the stronger is the tunnel rate. 
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FIG. 6. Scaling plot for the Laudau-Zener method showing 
the predicted field sweeping rate dependence for 1 T/s to 
1 mT/s. Each point indicates the magnetization after a field 
sweep over the M = ±10 resonance. The dotted lines are 
guides for the eyes. For comparison, there is also a relaxation 
curve without ac-field at H = 0. Inset: detail of the relaxation 
in the long time regime. 

We also compared the tunneling rates found by the 
Landau-Zener method with those found using a square- 
root decay method which was proposed by Pmkof 'ev and 
Stampo, and found again a good agreements. 

These measurements show that the Landau-Zener 
method is particularly adapted for molecular clusters be- 
cause it works even in the presence of dipolar and hyper- 
fine fields which spread the resonance transition provided 
that the field sweeping rate is not too smalltil. 
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7. Filed sweeping rate dependence of the effective 
tunnel splitting A measured by a Landau Zener method which 
works in the region of high sweeping rates where A is sweep- 
ing rate independent. The measured Landau-Zener tunneling 
probability P is indicated for several field sweeping rates. 



between the anisotropy hard axis and the transverse field 
(Fig. ||) show that for small ip angles the tunneling rate 
oscillates with a period between minima of ca. 0.41 T, 
whereas no oscillations showed up for large ip angles (see 
Fig. 2A in Ref.a). In the latter case, a much stronger in- 
crease of A with transverse field is observed. The trans- 
verse field dependence of the tunneling rate for different 
resonance conditions between the state M = —10 and 
(S—n) can be observed by sweeping the longitudinal field 
around H z = n x 0.22 T with n = 0, 1, 2, ..IB The cor- 
responding tunnel splittings A oscillate with almost the 
same period of ca. 0.4 T (Fig. ||). In addition, compar- 
ing quantum transitions between M = —S and (S — n), 
with n even or odd, revealed a parity (symmetry) effect 
which is analogous to the (Kramers) suppression of tun- 
neling predicted for half integer spinslprl. This behavior 
has been observed for n = to 4c3. A similar strong 
dependence on the azimuth angle ip was observed for all 
the resonances. 

In the frame of the simple giant spin model (Eq. 1), 
the period of oscillation (Eq. 2) is AH = 0.26 T for D = 
0.275 K and E = 0.046 K as in Ref.ll This is significantly 
smaller than the experimental value of ca. 0.4 T. In order 
to quantitatively reproduce the observed periodicity we 
have included forth order terms in the spin Hamiltonian 
(Eq.l) as recently employed in the swanulation of inelas- 
tic neutron scattering measurements^ and performed a 
diagonalization of the [21 x 21] matrix describing the S 
= 10 system. However, as the forth order terms are very 
small, only the term in C(S\ + St.), which is the most 
efficient in affecting the tunnel splitting A, has been con- 
sidered for the sake of simplicity. The calculated tunnel 
matrix elements for the states involved in the tunneling 
process at the resonances n = 0, 1, and 2 are reported 
in Fig. showing the oscillations as well as the parity 
effect for odd resonances. The period is reproduced by 
using D = 0.292 K E = 0.046 K as in RefE3 but with 
a different C value of -2.9 x 10~ 5 K. The calculated 
tunneling splitting, is however ca. 3 times smaller than 
the observed oneBJ. These small discrepancies are not 
surprising. In fact with the C parameter we take into 
account the effects of the neglected higher order terms 
in S x and S y of the spin hamiltonian which, even if very 
small, can give an important contribution to the period 
of oscillation and dramatically, affect A, as first pointed 
out by Prokof'ev and StampEj. Our choice of the forth 
order terms suppress the oscillations of A for |ffa;| > 1.4 
T which could not be studied in the current set-up. Fu- 
ture measurements should focus on the higher field region 
in order to find a better effective Hamiltonian. 



IV. OSCILLATIONS OF TUNNEL SPLITTING 

Studies of the tunnel splitting A, at the tunnel transi- 
tion between M = ±10, as a function of transverse fields 
applied at different angles ip, defined as the azimuth angle 
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FIG. 8. Measured tunnel splitting A as a function of trans- 
verse field for ~ 0°, and for quantum transition between 
M = — 10 and (S — n). Note the parity effect which is analo- 
gous ±0 the suppression of tunneling predicted for half integer 
spinstl. It should also be mentioned that internal dipolar and 
hyperfme fields hinder a quench of A which is predicted for 
an isolated spin. 
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FIG. 9. Calculated tunnel splitting A (Eq. 3) as a function 
of transverse field for quantum transition between M = — 10 
and (10 — n) at ip = 0°. Our choice of the forth order terms 
suppress the oscillations of A for \Hx\ > 1.4 T. Future mea- 
surements should focus on higher transverse fields. 

V. INTERMOLECULAR DIPOLE INTERACTION 

Fig. shows detailed measurement of the tunnel split- 
ting A around a topological quench for the quantum tran- 
sition between M — ±10, and M = — 10 and 9. Partic- 
ular effort were made to align well the transverse field 
in direction of the hard axis. The initial magnetizations 
5= M- m < M s were prepared by rapidly quenching the 
sample from 2 K in the present of an longitudinal applied 
field H z . The quench takes approximately one second 
and thus the sample does not have time to relax, either 
by thermal activation or by quantum transitions, so that 
the high temperature "thermal equilibrium" spin distri- 



bution is effectively frozen in. For H z > 1 T, one gets 
an almost saturated magnetization state. 

The measurements of A(M,- n ) show a strong depen- 
dence of the minimal tunnel spittings on the initial mag- 
netization (Fig. |Io| ) which demonstrates the transverse 
dipolar interaction between Fes molecular clusters being 
largest of A|i n = similar to the longitudinal dipolar 
interactior£3. 
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FIG. 10. Detailed measurement of the tunnel splitting A 
around a topological quench for the quantum transition be- 
tween M = —10 and (10 — n) at (fi = 0°. Note the strong 
dependence on the initial magnetization which demonstrates 
the dipolar interaction between Feg molecular clusters!^. 



VI. CONCLUSION 

Our measurement technique is opening a way of di- 
rectly measuring very small tunnel splittings of the order 
of 10~ 8 K not accessible by resonance techniques. We 
have found a very clear oscillation in the tunnel splittings 
A, which are direct evidence of the role of the topolog- 
ical spin phase in the spin dynamics of these molecules. 
They are also the first observation, to our knowledge, 
of an " Aharonov-Bohm" type of oscillation in a mag- 
netic system, analogous to the oscillations as a function 
of external flux in a SQUID ring. A great deal of infor- 
mation is contained in these oscillations, both about the 
form of the molecular spin Hamiltonian, and also about 
the dephasing effect of the environment. We expect that 
these oscillations should thus become a very useful tool 
for studying systems of quantum nanomagnets. 
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